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1 Introduction 

The simplest heat flow problem consists of a gas enclosed between two 
infinite, parallel plates kept at different temperatures. In the continuum 
limit, the Fourier law establishes a linear relation between the heat flux 
and the thermal gradient, i.e., q = — kVT, where k is the thermal conduc- 
tivity coefficient. The continuum description applies when A/L <C 1 and 
X/i <IC 1, where A is the mean free path, L is the separation between the 
plates, and £ ~ T| Vr|^^ is the characteristic length over which the temper- 
ature changes. Nevertheless, exact results from the Boltzmann equation for 
Maxwell molecules and from the BGK model for general interactions ^ ^ 
show that the Fourier law still holds for large gradients (A ^ £), provided 
that the Knudsen number Kn = X/L is small. 

In this paper we are interested in studying the influence on the heat flux 
of an external field g (e.g., gravity) normal to the plates. We will also 
assume that the flow velocity vanishes. This absence of convection is pos- 
sible if the Rayleigh number is smaller than its critical value (Ra < 1700) 
1^. This precludes the existence of the Rayleigh-Benard instability, that 
has been studied for dilute gases by other authors Q|. In addition to A, 
L, and £, the presence of gravity introduces a fourth characteristic length, 
namely h ^ ksT /mg, which represents the vertical distance over which the 
field produces a significant effect. In ordinary laboratory conditions, h is 
several orders of magnitude larger than A and £, so that the constitutive 
equations are not affected by the action of gravity. However, discrepancies 
with respect to the Navier-Stokes predictions can be expected if h is not ex- 
tremely large. According to a recent perturbation solution of the Boltzmann 
equation through second order in 5 [|| , one can estimate that the heat flux 
deviates from the Fourier law as much as 10% \i X ^ I ^ O.Ol/i. The aim of 
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this paper is to go beyond the second order in g by using the BGK model 
of the Boltzmann equation. Specifically, we will obtain the hydrodynamic 
profiles as well as the pressure tensor and the heat fiux through sixth order 
in gravity. 

2 Description of the problem 

Let us consider a dilute gas described by the BGK kinetic equation 

|/ + v.V/ + --|-/ = -;.(/-/i). (1) 
ot m ov 

Here, f(r,v,t) is the one-particle distribution function, F is an external 
force, m is the mass of a particle, i/ is a collision frequency, and Jl is the 
local equilibrium distribution function, that is characterized by the local 
density n, the local velocity u, and the local temperature T, defined as 

{n,nu,3nkBT} = J dv{l,v,m{v - uf}f, (2) 

where /cb is the Boltzmann constant. The collision frequency is proportional 
to the density and its dependence on the temperature models the infiuence 
of the interaction potential. For instance, for Maxwell molecules z/ oc n, 
while v cx nT^/'^ for hard spheres. 

The problem we want to investigate is that of a gas enclosed in a slab 
between two plates at different temperatures. We assume the existence of a 
stationary state with spatial variation along the normal direction z only and 
a constant external field F = —mgz along that direction. The constant g 
can be interpreted as the gravitational acceleration. In addition, we assume 
that there is no convection, i.e., u = 0. In order to ease the notation, it is 
convenient to introduce dimensionless quantities. To that end, we choose an 
arbitrary point zq in the bulk as the origin and take the quantities at that 
point (denoted by a subscript 0) as reference values. Therefore, we define 
T* = T/Tf), p* ^p/po, V* = v/uo, /* = {kBTo/pQ)v^f, g* = g/vQi^o, where 
p = nksT is the hydrostatic pressure and vq = (fcsTo/m)^/^ is a thermal 
velocity. In these units, g* = Xo/ho, where we define the mean free path as 
•^0 = i^o/fo- In the case of the spatial variable z, it is convenient to rescale 
it in a nonlinear way that takes into account the local dependence of the 
collision frequency. Consequently, we define 

s — / dz'v{z'). (3) 

J Zq 

Thus, the stationary BGK equation reads 

l + <a.-.g*^i?„)r (4) 
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where ds = d/ds and Dy = d/dv^. Furthermore, for the sake of concrete- 
ness, we have restricted ourselves to the case of Maxwell molecules (i.e., 
V oc p/T). In the geometry of the problem, the relevant velocity moments 
are defined as 

M„^= [ d^r^v^^'^vfr. (5) 



In particular, the corresponding moments at local equilibrium are 

M^^ = (^^±A±l)!lp*r— i+^'Z^ (6) 

for (3 even, being zero otherwise. Conservation of momentum and energy 
implies that 9sMo2 — —g* and dgMn = 0. 

In the absence of gravitation {g = 0), Eq. has an exact solution |^ 
characterized by a constant pressure, p* = 1, and a "linear" temperature 
profile, T* = 1 + es, that applies to arbitrary values of the reduced thermal 
gradient e — Xo/£o. The velocity moments are polynomials in both s and e. 
Their explicit expression for + 2(a — 1) > is Q 

(a-l + ^)!(/3 + r + l) 

(r+/3)even 

(7) 

In particular, Mn = — 5e, which means that the Fourier law holds even for 
large thermal gradients. 

The motivation of this paper is to analyze the influence of gravitation on 
the profiles and transport properties of the above steady Fourier flow. How- 
ever, the presence of the operator Dy in Eq. (^) complicates the problem 
significantly. A convenient strategy is to take the pure steady Fourier flow 
corresponding to a value of e equal to the actual thermal gradient at the 
point s = as a reference state. Consequently, we will carry out a pertur- 
bation expansion in powers of g: 

= + M^a* + M^g*' + • • • , (9) 
p* =p^"'> +p^'^g* +p^^'>g*^ + ■■■, (10) 

T* = r(0) + + T^2)g*2 ^ , , , ^ ^^^^ 

where M^^^ is given by Eq. (^, p'^^ 1, and T*^"^ — 1 + es. By definition, 
p('=)(0) = T('=)(0) = dsT'-''^^^^ for fc > 1. It must be emphasized that 
the terms of order g*'^ are nonlinear functions of e since no restriction to 
the order on e exists. 
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3 Perturbation expansion 

In this section we obtain the hydrodynamic profiles p'-*^'' and T*-*^-*, the mo- 

(k) (k) 

mentum flux Mq2 , and the heat flux through order k — 6. Insertion 

of Eq. (|) into Eq. (|) yields 

OO 

j=o 

This is a formal solution, since f'^^ is a functional of f^''^ through its de- 
pendence on the pressure and temperature. Taking moments in Eq. (]l^), 
one has 

j=l j=0 i=0 ^ ^ ^ 

X {^^Mtr,i]+i + (/? + ^■)<^":7-i) ■ (13) 

The fact that / and /l have the same hydrodynamic moments leads to the 
consistency conditions 

AAf^o^ = AM^'l^ = Ail/^o' = 0, (14) 

for any k. In order to convert Eq. (|l^) into an explicit equation that can be 
solved recursively, we need to know the spatial dependence of p'^'^^ and T'^'^^ . 
It turns out that p^*"'' is a polynomial of degree fc — 2 in s (except p''^\ that 
is linear in s), while T*^'^-' is a polynomial of degree fc + 1, the coefficients 
being nonlinear functions of the reduced thermal gradient e. Thus, at a 
given order fc, insertion of these poly nomials into Eq. (^|) and application 
of the consistency requirements ( |l4[) allow one to determine the unknown 
coefficients and the problem can be recursively solved. Notice that, seen 
as functions of the actual space variable z, p^'''> and T^''^ are much more 
complicated than just polynomials. In fact, in the absence of gravitational 
force, the relationship between s and z is nonlinear: z = Xo{s + ^es^). In 
general, such a relationship can be obtained inverting Eq. (jsj) as 

z = ZQ + XoJ^ d^'^;^- (15) 

The above scheme is straightforward but tedious to carry out. Since all the 
manipulations are algebraic, they render themselves to the use of symbolic 
programming languages. In this paper, we have evaluated the perturbation 
expansion through sixth order in the field. Since the expressions of the 
hydrodynamic profiles become progressively longer, here we only give the 
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explicit results through order k — A: 



, , 276 2 ,3 1 
1-S5 e sg --s 



— e (112 9736^ + 30) + 588e^s 
5 



(16) 



= l + es + -es^g* 

346 1 , 
15 4 



66 
5" 



, *2 2 



(66246^ + 5) 



5" 



i^e(50 765 962e2 + 31445) 



2 971 

(399 621e2 + 200) s + es^ - 

lb '6 



g** + 0{g*'). (17) 



Once the hydrodynamic profiles are known, Eq. (|13[ ) can be used to obtain 
all the velocity moments at a given order. The most relevant moment is Mn, 
which is related to the heat flux, = {pqVq/2)Mii. Another important 
quantity is Afo2 = Pzz/Poi which measures the anisotropy of the pressure 
tensor P. As said above, in the absence of gravitation the Fourier law applies 
exactly, i.e. gi"-* — —ndT/dz — — (5po'i'o/2)e, and the pressure tensor is 
isotropic, i.e. pj? = po- In order to characterize the deviations from these 
Navier-Stokes predictions due to gravity, we introduce the following reduced 
quantities: 



A(e,g*) = ^ = l + 5]AW(6)g*\ 



(18) 



k=l 



Pa 



1 + J2^^'\e}g*\ 



(19) 



fc=i 



The results show that A^^*^) and 7'^'"'-' are polynomials in e of degree k and a 
defined parity. 



A«(e) = ^Af^.^ 7W(,) = ^^«,^ 



(20) 



To second order in g* we get A^^' = ^, A^,^^ = ^, A^^^ 



5 ' '"-2 



47 968 (1) _ Q 
25 ' '1 ' 



5 ' '"-0 

and 72^"* = ^ The coefficients A^*^'' and 7^'''' for 3 < fc < 6 are given in 



Tables |l| and^, respectively. 



4 Discussion 

The numerical coefficients appearing in Eqs. (H) and (0), as well as in 
Tables |l| and || clearly indicate that the expansion in powers of g*, Eq. (^), 
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Table 1: Coefficients lS.^f^ for 3 < fc < 6 
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5 25 
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^ 1579 897 956 43 675 137132 624 

125 625 
20 495 361514 332 



625 



449 964 261460 472 196 
ai2S 



Table 2: Coefficients 7^ ^ for 3 < /c < 6 



is only asymptotic. This does not pose a serious problem, at least from a 
practical point of view, except for values of g* (say, g* > 10~^) that cor- 
respond to gravitational fields unrealistically large. As a consequence, only 
the first few terms in the expansion are useful for small values of g* . In the 
subsequent analysis, terms of order g*^ and higher will not be considered. 
This also allows us to make a closer comparison with results derived from 
the Boltzmann equation for Maxwell molecules Notwithstanding, the 
knowledge of the remaining terms might be useful to attempt to resum the 
infinite series and get the transport properties for arbitrary values of g*. 
The solution to the Boltzmann equation |^ through order g*^ exhibits the 
same structure as in Eqs. (|l^)-(p^), so that only the numerical coefficients 
differ. As a matter of fact, the coefficient ^ appearing in Eq. ( p7| ) is replaced 

by in addition, the Boltzmann solution yields A^^-* = Ag^^ — 

A^^^ = 503.7, 7^' = and j^^'' = The differences indicate that the in- 
fluence of gravity is stronger in the BGK description than in the Boltzmann 
one, especially in the case of the pressure anisotropy. In order to carry out a 
more detailed comparison, it is convenient to construct Fade approximants 
for the generalized thermal conductivity coefficient A. More specifically, we 
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consider the approximants 



A[i,i](e,5*) 



A(i)(e) + 


A(i)'(e)-A(2)(e)' 


9* 


A(i 


(e)~A(2)(e)5* 



A 



[0,2] 



(e,/) = {l - A^'\e)g* + [A^'^\e) A^'\e)] g*'} 



(21) 



(22) 



In the case of the BGK equation, these two approximants differ less than 



1.075 




Figure 1 : Plot of the (reduced) nonlinear thermal conductivity A versus the 
thermal gradient e ai g* = 10~^, as obtained from the BGK model (solid 
line) and the Boltzmann equation (dashed line). 



2% if g* < 10~^ and \e\g* < 3 x 10^''. The differences are smaller in the case 
of the Boltzmann equation. In this range of values for e and g*, a reliable 
approximation is A ~ i (^[1,1] + A[o 2])- Figure [l] shows this approximation 
for A in the interval —3 < e < 3 at 5* = 10"'^, as given by the BGK and 
Boltzmann equations. We observe that the heat flux increases with respect 
to its Navier-Stokes value when one heats from above (e > 0), while the 
opposite happens when one heats from below (e < 0). This effect is not 
symmetric, since it is more significant if e > than if e < 0. As said 
above. Fig. |^ also shows that the influence of gravity is less important in 
the Boltzmann description. 

In summary, we have solved the BGK model for a gas simultaneously sub- 
jected to a thermal gradient and a parallel gravity field of magnitude g, in 
the absence of convection. The solution has been obtained by a perturbation 
expansion in powers of gravity, the reference state being the non-equilibrium 
pure Fourier flow with arbitrarily large thermal gradients. We have explic- 
itly obtained the solution through sixth order in the field. The results clearly 
indicate that the expansion is not convergent, although it seems to be at 
least asymptotic. The work reported in this paper extends previous results 
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derived from the Boltzmann equation to second order in the field Q. The 
similarity in the structure of the coefficients appearing in both descriptions 
suggests that the expansion obtained from the Boltzmann equation is also 
asymptotic. Nevertheless, given that at practical level, the values of g are 
small, the usefulness of the expansion is restricted to the first few terms. 
The main results concerning the transport of momentum and energy are that 
the external field induces (i) anisotropy in the pressure tensor, {Pzz —p)/p — 
^e^g*^, and (ii) deviations from the Fourier law, q^jiiz^ — 1 ~ ^tg* . 
While the first effect is of second order, the correction to the heat flux is 
of first order, so that it depends on the sign of the thermal gradient. As 
a consequence, the heat transport is inhibited when the gas is heated from 
below (e < 0), while the opposite happens when the gas is heated from 
above (e > 0). 
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